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Abstract 

Ginzburg-Landau model with two order parameters appears in many condensed- 
matter problems. However, even for scalar order parameters, the most general U{\)- 
symmetric Landau potential with all quadratic and quartic terms contains 13 indepen- 
dent coefficients and cannot be minimized with straightforward algebra. Here, we develop 
a geometric approach that circumvents this computational difficulty and allows one to 
study properties of the model without knowing the exact position of the minimum. In 
particular, we find the number of minima of the potential, classify explicit symmetries 
possible in this model, establish conditions when and how these symmetries are sponta- 
neously broken, and explicitly describe the phase diagram. 

1 Introduction 

The Landau theory [TJ [2] offers a remarkably economic description of phase transitions as- 
sociated with symmetry breaking. This breaking is described by an order parameter ip: the 
high symmetry phase corresponds to ip = 0, while the low symmetry phase is described by 
ip 7^ 0. Very often, the order parameter can be directly related to physically observable quan- 
tities, such as, for example, distortion of the crystal lattice or spontaneous magnetization, The 
local version of the Landau theory with a coordinate-dependent order parameter, known as 
Ginzburg-Landau (GL) theory, is the basis of the phenomenological theory of superconduc- 
tivity, [3]. For a variety of applications of the Landau theory to various condensed-matter 
problems, see e.g. [I] and references therein. 

In order to find if a given system is in its high or low symmetry phase, one constructs a 
Landau potential that depends on the order parameter, and then finds its minimum. For a 
complex order parameter, its classic form [TJ, [2] is 

VW = -a\^+ b -W + o(\^). (1) 
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Near the phase transition the higher order terms o(|^| 4 ) are usually assumed to be negligible. 
The values of the coefficients a and b and their dependence on temperature, pressure, etc. can 
be either calculated from a microscopic theory, if it is available, or considered as free parameters 
in a phenomenological approach. The phase transition associated with the symmetry breaking 
takes place when an initially negative a becomes positive, and the minimum of the potential 
(TjQ) shifts from zero to 

W = y|e fa , (2) 

with an arbitrary phase a. 

Many systems are known in which two competing order parameters (OP) coexist. Among 
them are the general 0(m) © 0(n)-symmetric models, [5], the models with two interacting N- 
vector OPs with 0(N) symmetry, [6]; spin-density- waves in cuprates, j7]; competition between 
antiferromagnetism and superconductivity, [8] ; 4 He with its interplay of crystalline and super- 
fluid ordering, [9]; multicomponent, pU], non-conventional two-dimensional, [llj, spin-triplet 
p-wave, p2], and two-gap, [131 E]> superconductivity, with its application to magnetism in 
neutron stars, [15]; two-band superfluidity, |16| ; and even mechanisms of electroweak sym- 
metry breaking beyond the Standard Model such as the two-Higgs-doublet model (2HDM), 

m 

To describe such a situation within GL theory, one constructs a Landau potential similar 
to ([T]), which depends on two order parameters, fa and fa. With scalar order parameters, it 
can be written generically as 

V(fa,fa) = -a^tfa) + -b ijk i(ip*ipj)(iplfa) , k, I = 1, 2 . (3) 

The coefficients of this potential can be considered independent although in each particular 
application they might obey specific relations. One thus arrives at the general two-order- 
parameter (20P) GL model with quadratic and quartic terms. 

Once Landau potential ([3]) is written, the next step is to find its minimum, i.e. to solve 
the static homogeneous Ginzburg-Landau equations. A rather surprising fact is that these 
equations cannot be solved with straightforward algebra. Differentiating the Landau potential 
with respect to fa leads to a system of coupled algebraic equations of third order, whose total 
degree of algebraic complexity is six, which makes it impossible to solve in the general case. 

In this paper we argue that despite this computational problem, there remains something 
that one can learn about the most general two-order-parameter model in the mean-field ap- 
proximation: its phase diagram. As we will show below, it is possible to classify all the phases 
according to the symmetries of the model and properties of the ground state. 

This idea is not new. In fact, there exists an extensive literature dating back to 1980 's on 
minimization of G-invariant potentials with several OPs (G being a group of transformations 
of OPs), see e.g. [U [18], [19] and references therein. These works exploit the fact that the 
problem becomes simpler when reformulated in the orbit space instead of the space of order 
parameters themselves, [20J. This orbit space is naturally sliced into several strata, which are 
linked to the allowed phases of the model. To describe them, one constructs the ring of G- 
invariant polynomials of the order parameters and finds the Minimal Integrity Basis (MIB) of 
this ring. Different strata (i.e. different phases) correspond to some particular relations among 
the MIB polynomials. This general method has been applied to classification of the phases 
of several relevant physical systems, for example, to p-wave superfluidity of 3 He, [21], D-w&ve 
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condensates, [22], 2HDM, [231 [21], and to general Landau models with multicomponent and 
even multidimensional order parameters, [25, 26]. 

In the light of this activity, it is somewhat surprising that the most general GL model 
with two complex order parameters and with the most general potential in the form of ([3]) 
has never been studied in complete detail. Here we fill this gap showing that in this case the 
analysis can be pushed much farther than in the general situation, with important physical 
consequences. 

The approach presented here is based on the reparametrization symmetry of the model, 
which allows one to establish the Minkowski-space structure of the orbit space. The minimiza- 
tion problem admits a transparent geometric interpretation, which leads to several theorems 
concerning the properties of the global minimum. Specific application of this approach to the 
2HDM was given in [28]. [29] . Here we analyze the case of two local order parameters ipi(r) and 
?/>2(^*) in the general context, which can be relevant also for many condensed-matter problems. 

1.1 Geometric approach vs. Minimal Integrity Basis method 

Let us stress from the very start the essential differences between the geometric analysis of 
the present paper and the standard approach based on the Minimal Integrity Basis (MIB) 
technique. 

The first difference lies in the scope of these two approaches. The MIB leads to interesting 
results in the cases when the potential is invariant under a non-trivial group G of transforma- 
tion of the order parameters. The larger G, the richer is the spectrum of possible patterns of 
its spontaneous violation. In particular, MIB methods has nothing to say if G is the trivial 
group. 

In a typical situation one takes a highly symmetric G-invariant potential constructed from 
powers of several multidimensional order parameters up to a certain degree, builds various in- 
variants, finds the ones that form MIB, and classifies the possible phases according to relations 
among these invariants. This approach is rather general in a sense that it can be applied, in 
principle, to any number and any dimensions of order parameters. However, because of G- 
invariance, the potentials usually contain very few terms. 

The geometric approach presented here is limited to the particular case of two complex 
scalar or vector order parameters, and to the fourth-degree potentials. However, within these 
restrictions, we manage to work out the most general model with all possible types of the 
OP interactions. The only symmetry that we impose is the £7(l)-symmetry of the free energy 
density, which is a reasonable choice from the physical point of view. In this aspect, our 
analysis is more general than the MIB approach: we just take two OPs, construct the free 
energy density in its full complexity, and study everything that can ever happen in this model. 
That is, we analyze all possible symmetry groups G and all possible patterns of symmetry 
breaking. 

The second difference concerns the procedures and the results of these two approaches. In 
a situation when several phases are possible, one wants to know which phase corresponds to 
the ground state of the model (i.e. which phase is stable). In the usual MIB method one can 
do nothing but explicitly solve the algebraic equations and check the minimum conditions. 
This can be done only if the equations are simple enough, which in turns happen when the 
free energy density is simple. Thus, only sufficiently symmetric potentials are fully tractable 
with the MIB method. Examples cited in [2_H |2B] are precisely of this type. 
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In terminology suggested by [27], one should distinguish between the angular problem 
(classifying all the allowed phases) and the radial problem (actually finding the position of the 
absolute minimum of a given potential). MIB methods allow one to solve the angular, but not 
the radial problem. 

In the case of the most general 20P GL model, with its large number of free parameters, 
this algebra cannot be worked out explicitly. One ends up with a general algebraic equation 
of sixth order, which one cannot solve analytically. Thus, one is unable to solve the problem 
of minimization of the potential using only the MIB formalism. 

The strongest point of the present geometric approach is that we avoid solving these equa- 
tions and nevertheless we rigorously prove several statements about the ground state of the 
model. In other words, we study the properties of the absolute minimum without solving the 
radial problem. This is especially useful for the case of the smallest possible group G, for 
which the MIB technique becomes redundant. Thus, the geometric approach presented here 
is neither a particular case nor an improvement of the MIB method, but is complementary to it. 

The structure of the paper is the following. In Section [2] we introduce the formalism and 
derive a very compact expression for the free energy functional. The extrema of the Landau 
potential cannot be found with straightforward algebra, so in Section [3] we develop geometric 
tools that allow us to find the number of extrema and minima of the potential. Section H] 
is devoted to the special case of a potential stable in a weak sense. Then, in Section |5] we 
give full classification of explicit symmetries of the model and derive conditions when and 
how these symmetries are spontaneously broken. All this allows us to describe in Section [6] 
the phase diagrams of the model, listing the phases according to the number of minima and 
symmetries. Here, we also discuss phase transitions and argue that critical properties, too, can 
be calculated in geometric terms. Section [7J contains analysis of several simple cases, which 
provide illustration of the general approach. In the short Section [H] we outline conditions when 
solitons appear in this model. In Section we outline characteristic features of the general 
GL model with two complex iV-vector order parameters, and finally in Section [10] we draw 
our conclusions. Appendices provide some mathematical details and derivations. 

2 Formalism 

In the main part of the paper we will assume that ipi{r) are just complex numbers; modifi- 
cations in the case of more complicated OPs will be discussed in Section [9j Throughout the 
paper we also assume that the absolute values of 1^1 are not bounded from above. 

Let us consider the free-energy density in the most general globally [/(l)-invariant 20P 
GL model containing all possible quadratic and quartic terms in the potential: 

F = K + V 2 + V A . (4) 

It is a sum of the gradient term K, 

K = k x \D^\ 2 + k 2 \D^ 2 \ 2 + k 3 {DiM*{D^) + 4(D^ 2 Y(D^) , (5) 
where D is either V or the covariant derivative, and the Landau potential 

V 2 = -ai|^i| 2 - a 2 \^ 2 \ 2 - 03(^1^2) - aKip^i) 5 (6) 
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V, 



\H 4 +-^\H 4 +bM 2 \H 2 + 



(V^) + C.C. 



Free energy density (TjJ contains 4 + 4 + 9 = 17 free parameters: real k±, k 2 , a%, a 2 , b\, b 2 , &3 
and complex k 3 , 03, 64, 65, b & . 

By construction, the free energy remains invariant under the U(l) group of simultaneous 
multiplication of ip± and ip2 by the same global phase factor. We do not consider terms that 
violate this symmetry, such as ipf + ijpi)* ■ 

Note that potential contains quartic terms such as | ^1 1 2 C V^i ^2 ) that mix ipi and ip2, 
which are usually absent in many particular applications of the 20P GL model. However, in 
certain cases such terms appear, as it happens in the dirty limit of a two-gap superconductor, 
see e.g. [14]. 

We stress that in our approach it is essential that we include all possible terms from the 
very beginning. 



2.1 Reparametrization symmetry 

From the physical point of view, the order parameters tpi and ip 2 can be of the same (as in two- 
gap superconductors) or of different nature (as in the case of superfluid/crystalline ordering 
interplay). However, one can always make OPs dimensionless, and once the free energy density 
(T4J) is constructed and the problem of its minimization is posed, the physical nature of the 
OPs becomes irrelevant. 

One can then view OPs ip\ and ip 2 as components of a single complex 2-vector $: 

-(*)■ 

and consider transformations that mix ipi and ip 2 . These are assumed to be local transforma- 
tions, i.e. they mix ipi{r) taken at the same point r. 

We start with the observation that the most general free energy density (jl]) retains its 
generic form under any regular linear transformation between ipi and ip 2 . In other words, 
transformation 

$ _> $' = T ■ $ , with any T e GL{2, C) , (7) 
again leads to (jlj) but with reparametrized coefficients: 

{Ki, a t , h} -> {k-, a-, &•} = r(«i, a i? k) . (8) 

The explicit link between T and r will be given below. 

Since any T 6 GL(2,C) is invertible, so is r. Therefore, if (j7]) is accompanied by the 
transform r _1 of the coefficients, then one arrives at exactly the same expression for the free 
energy as before. 

If one considers the free energy only, then the physical observables, such as the depth of 
the Landau potential at the minimum and the eigenvalues of the second derivative matrix 
of the potential (the hessian) can be expressed in terms of the coefficients {Ki, ai, b{\ only. 
Therefore, the models ($, {Ki, a iy bi}) and (<&', a[, b'j}) have the same sets of observables. 
In other words, reparametrization transformations do not change the physical content of a 
given model; they only affect the way we look at it. Thus, we have a reparametrization 
freedom in this problem, with the reparametrization group GL(2, C). 
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The general linear group GL(2, C) is an 8- dimensional Lie group. It can be written as 



GL(2,C) = C* x SL(2,C), 



(9) 



where C* is the group of all multiplications of $ by a non-zero complex number. Due to the 
[/(l)-invariance of the free energy, multiplication of $ by an overall phase factor induces the 
identity transformation of the coefficients, while the 7-dimensional factorgroup GL(2, C)/U(l) 
induces non-trivial transformations r. Thus, the 1 7-dimensional space of coefficients (i.e. the 
space of all possible 20P GL models) becomes sliced into 7-dimensional regions of essentially 
identical models linked by all possible r. The space of distinct physical situations is described 
by the corresponding 10-dimensional factorspace. 

2.2 Orbit space 

Let us now introduce the four-vector r M = (ro, r^) = ($V M $) with components 



Here, index \i = 0, 1,2,3 refers to the components in the internal space and has no relation 
to the space-time. Multiplying ipi by a common phase factor does not change r M , so each r M 
uniquely parametrizes a single [/(l)-orbit in the ^-space. The ?7(l)-invariant free energy (Tj0) 
can be also defined in this 1 + 3-dimensional orbit space. 

The SL(2, C) C GL(2, C) group of transformations of $ induces the proper Lorentz group 
50(1, 3) of transformations of r M . This group includes 3D rotations of the vector as well as 
"boosts" that mix r and rj, so the orbit space gets naturally equipped with the Minkowski 
space structure with metric diag(l, —1, —1, —1). We stress again that the words "Minkowski 
space" and "Lorentz group" always refer to the internal space, not to the usual space-time. 

Since the order parameters ipi an d V>2 are just complex numbers, direct calculation shows 
that 



Then, since r > and since the values of r M are not restricted from above, the orbit space of 
the 20P GL model is given by the forward lightcone LC + in the Minkowski space. As should 
be expected, the reparametrization group in the orbit space, SO(l,3), leaves the orbit space 
invariant. 

Analogously to r M , one can also introduce 




2Re (^2) 
21m (^2) 

\H 2 -\H' 




(10) 




(ii) 



(12) 



Obviously, the reparametrization transformation laws of are the same as for r M . 
All this allows us to rewrite the free energy (j3J) in a very compact form: 



(13) 
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with 



B 



fJtU 



- («! + k 2 , -2Re k 3 , 21m k 3 , —k x + k 2 ) 



- (ai + a 2 , — 2Re a 3 , 21m a 3 , — ai + a 2 ) 



V 



61+62 I l. 
—3— + 3 

-Re (65 + fe 6 

Im (65 + fe 6 ) 
61—62 



-Re (65 + fee) 
Re fe 4 
— Im fe 4 
Re (65 - fe 6 ) 



Im (65 + fe 6 ) 
— Im fe 4 
—Re fe 4 

-Im (65 - fe 6 



-62 



Re (fe 5 - fo 6 ) 
-Im (65 - fe 6 ) 

—5 °3 / 



(14) 



Note that due to l ull) , definition of the matrix has one degree of freedom, since B^ v and 
B^ v = B^ v — Cg^ u with any C are equivalent. 

The quantities K^, A 11 and B^ u transform as four- vectors and a four-tensor, respectively. 
This provides the explicit link between transformations T and r mentioned in Sect. 12.11 For 
convenience, we collect in Appendix |A] some basic facts concerning the manipulation of B^ u . 



2.3 Properties of the coefficients 

General physical requirements place restrictions on possible and B^ v . 

First, the requirement that very large wavevector oscillations must increase not decrease 
the free energy implies that lies inside the future lightcone: K$ > 0, K^K^ > 0. This 
condition remains true under an arbitrary 5*0(1,3) transformation. 

Second, we require that the potential is bounded from below in the entire ^-space. Since 
the potential is build of quadratic and quartic terms, V — V 2 + V4, this can be achieved in 
two cases (here we coin the terminology of [30], where the stability of the Higgs potential in 
2HDM was analyzed): 

• the potential is stable in a strong sense, if V4 increases in all directions in the ^j-space; 

• the potential is stable in a weak sense, if V4 is non-decreasing in all directions in the 
Vvspace, and V 2 increases along the flat directions of V4. 

Let us focus on the case of the potentials stable in a strong sense; the case of the potential 
stable in a weak sense will be considered in Section HI The requirement that V4 is positive 
definite in the entire ^j-space means that the quadratic form B^r^r" is positive-definite on 
the future lightcone LC + . In Appendix [B] we prove that this is equivalent to the statement 
that B^ is diagonalizable by an SO(l,3) transformation and after diagonalization it takes 
form 



B 



fJ,U 



( B 
















-B x 
















-B 2 







^ 








-B 3 


J 



with B >Bi,B 2 ,B 3 . (15) 



We will refer to B as the "timelike" eigenvalue of B^ v and B iy i = 1,2, 3, as its "spacelike" 
eigenvalues. The sign minus in front of the spacelike eigenvalues is the result of the Minkowski- 
space metric, see Appendix [A] The degree of freedom in the definition of B^ v amounts to 
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shifting all the eigenvalues by the same constant and does not affect the inequalities (IT5|) . 
However it can be used to manipulate the signs of the eigenvalues. 

Finding the eigenvalues of B^ v explicitly in terms of 6j requires solution of a fourth-order 
characteristic equation, which constitutes one of the computational difficulties of straightfor- 
ward algebra. We reiterate that in our analysis we never use these explicit expressions. Our 
analysis relies only on the fact that the eigenvalues are real and satisfy ( TT5l ). 

3 Minima of the Landau potential 

Having introduced the formalism that allows us treat the most general 20P GL model, let us 
proceed to the task of minimization of the free energy functional. We do not consider here 
the effects of non-trivial boundary conditions, so we are looking for homogeneous solutions 
ipi{r) = (ipi) that minimize the Landau potential ([6]). 

As mentioned in the Introduction, straightforward algebra is of little help for the minimiza- 
tion problem, since the resulting system of coupled equations dV/ dipi = cannot be solved in 
the general case. However one can still learn much about the ground state of the general 20P 
GL model without finding its location explicitly. In this paper we will provide, in particular, 
answers to the following questions: 

• How many extrema does the potential with given parameters have? How many of them 
are minima? 

• Can the global minimum be degenerate and when does it happen? 

• When does the global minimum spontaneously break an explicit symmetry of the po- 
tential? 

• What is the phase diagram of the model? What phase transitions can take place during 
continuous change of the coefficients of the model? 

3.1 Number of extrema 

Let us start with the number of extrema of a generic Landau potential. In order to find 
an extremum of V lying on the future lightcone LC + , one can use the standard Lagrange 
multiplier method. In this case one needs to introduce only one Lagrange multiplier A, which 
leads to the following system 

/ B, u (r»)-\(r,) = A,, 

\(r ft )(r lt ) = 0. [ ' 

Here, (r^) labels the position of an extremum. To avoid cumbersome notation, we omit (. . .) 
in this subsection. 

To establish how many solutions system ( lT6i) has, consider the .B^-diagonal frame (we 
remind that for a potential stable in a strong sense such a frame always exists). Then the first 
line in (fT6|) takes form 

(B - A)r = A , {Bi - X)n = A, . (17) 
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Rewriting = r nj, where is a unit 3-vector, and eliminating A, one obtains 



[A Q - (B - B t )r } n i = A i . (18) 

These three equations are coupled via the condition |n| = 1. Consider the l.h.s. of (1181) at 
fixed r and all the unit vectors n.j. It parametrizes an ellipsoid with semiaxes 

A -(5o-£i)r , A -(B -B 2 )r , A - (B - B 3 )r . (19) 

Now imagine how this ellipsoid changes if r increases from zero to infinity. Let us for simplicity 
assume that the eigenvalues of B^ v are distinct and Bi < B 2 < B 3 . 

Assume first that A > 0. Then, at r = 0, Eq. ( fl8l) parametrizes a sphere with radius A . 
As ro increases, it turns into a continuously shrinking ellipsoid with semiaxes ( fl9l) . At 

(i) _ A 
r = rX' 



B — £?x 

this ellipsoid collapses to the interior of a planar ellipse with semiaxes 

B 2 — B\ B% — B\ 

-DO — -t>0 — -Dl 

orthogonal to the first axis. As r increases further, this ellipse returns to an ellipsoid with 
two shrinking and one growing semiaxes, and at 

(2) A) 

r ° = r ° = R R" 

-DO — -D2 

is collapses again to a flat ellipse with semiaxes 

I B\ — B2 1 . -B3 — B2 

0_ R R ' Aq ~R W 1 

-DO — -D2 -DO — -D2 

(3) 

orthogonal to the second axis. Further on, at r = r Q it collapses to an ellipse orthogonal to 
the third axis, and for even larger values of r this ellipsoid increases infinitely. 

For each r interval, the ellipsoid sweeps a certain region in the three-dimensional space. 

• During the first stage, < r < r^, it sweeps the interior of the sphere of radius A , 
passing through each point exactly once. 

• during the second stage, < r < it sweeps a certain region, bounded by the 
caustic surface shown in Fig. [TJ left. It can be shown that each point inside this region 
is swept exactly twice. 

(2) (3) 

• during the third stage, r ; < r < r , it sweeps twice a similar caustic region, but 
oriented differently; 

(3) 

• Finally, during the fourth stage, r > r , it sweeps once the entire 3D-space. 
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Figure 1: (Color online) Left: The envelope of ellipsoids for r$ < r < Tq. Right: (A\, As- 
sertion of the caustic surfaces in the Aj-space and of the sphere with radius Aq. The number 
of solutions of Eq. ffT8l) is indicated for each region. 

Note that the appearance of caustic regions in the potential extremization problem is natural, 
since this problem is known to exhibit some catastrophe theory phenomena, see e.g. 0j. 

Returning to the system (fl8l) . which equates the l.h.s. to the 3- vector Ai, one sees that in 
order get the number of solutions of (ITS]) without finding them explicitly, one simply has to 
check whether Ai falls inside these regions. Fig. [TJ right, illustrates this statement. It shows 
regions with different numbers of extrema on the {A\, A3) plane for A 2 = and some Aq > 0. 

For the second possibility, A < 0, the situation is much simpler. At r = we again start 
with the sphere of radius \A \. As r grows, it turns into an ellipsoid with growing semiaxes 
(which is due to B — Bi > 0), and it sweeps once the entire space outside the sphere. 

The size of the above 3D regions is proportional to A . Therefore, in the 1 + 3-dimensional 
space of four- vectors A^, they define the corresponding conical regions starting from the origin. 
Therefore, the number of extrema of the potential depends on where the four-vector A 11 lies: 

• If A^ lies inside the past lightcone LC~ (i.e. Aq < and \A\ < \Aq\), then system ( fl6i) 
has no solution. In this case the quadratic term of the potential, —A^r 11 , increases in all 
directions in the ^-space. The only extremum of the potential is the global minimum 
at the origin, which corresponds to the high-symmetry ground state of the model. 

• If A 11 lies outside LC~ , then at least one non-trivial solution exists. If Aq < in the 
.B^-diagonal basis, i.e. A^ still lies in the lower hemispace, then this solution is unique 
and is the global minimum of the potential. 

• If A M lies inside LC + , then at least two non-trivial extrema exist. 

• If A^ lies inside one or both caustic cones defined above, then two additional extrema 
per cone appear. 

In total, there can be up to six non-trivial extrema of the potential in the orbit space. This 
result was also found independently in [30] with a more traditional analysis of the Higgs 
potential of 2HDM. The largest number of extrema is realized in situations when Aq > and 
Ai is sufficiently small, so that A^ lies inside both caustic cones. 
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Special care must be taken when r of an extremum is exactly equal to one of the values 
70 . Then the sequence of intersections of the ellipsoid with a given point changes, but 
the overall counting rules given above remain the same. As we will see later, this situation 
corresponds to spontaneous violation of a discrete symmetry. 



3.2 Number of local minima 

In general, the above construction cannot distinguish a local minimum from a saddle point or 
a maximum, so other methods must be used to establish the number of local minima. 

First of all, let us note that potential ([6]) with restrictions (1151) cannot have non-trivial 
maxima [20, 28J. This can be easily seen by drawing any ray in the ^-space from the origin 
and observing that the potential along this ray can be written asa;|'0| 2 +/3|?/>| 4 with (3 > 0. This 
function can never have a non-trivial maximum. Thus, the problem reduces to distinguishing 
minima from saddle points (in the orbit space). 

Take a generic extremum of the potential in the ipi, ip2 space, and calculate the second 
derivative matrix of the potential (the hessian) at this point 

= irrTrr ■ (20) 



Here, Q are the four real degrees of freedom, real and imaginary parts of ipi and ^2 : 

ipx = (pi + i(f) 2 ; ip 2 = 03 + #4 • (21) 

We will refer to the eigenvalues of this matrix as "eigenfrequencies" , uj 2 a . Due to the £7(1)- 
invariance of the potential, it always has one flat direction with zero eigenfrequency (one 
Goldstone mode), while among the other three, there is at least one positive eigenfrequency. 

Let us call the signs of these three eigenfrequencies (i.e. + + +, +H — , or H ) the signature 

of the hessian. 

Among the four degrees of freedom in the i/'i-space, three correspond to variations in the 
orbit space, i.e. to shifts of the point r 1 ' on LC + away from the extremum. If the extremum 
is not at the origin, then these shifts are linear functions of the shifts in the ^-space, and the 
Jacobian corresponding to this transformation is regular. Indeed, with the notation ([2"Tj) . one 

gets: 

/ 01 02 03 04 \ 

(22) 





/ 


01 02 03 


04 \ 


1 dr^ 




03 04 01 


02 


2<90 a 




04 -03 -02 


01 




V 


01 02 -03 


-04 / 



If $1 is non-zero, then this matrix has one and only one zero eigenvalue, with the corresponding 
right eigenvector (—02, 0i, —04, 03) being the Goldstone mode. This can be seen most easily 
in the frame where oc (1,0,0,1), implying 3 = 4 = (obviously, such a frame always 
exists for any r^). The signature of the hessian, therefore, is the same in the ^-space and in 
the orbit space. 

For an extremum to be minimum, its signature must be + + +. However, since the 
explicit expressions for (r^) cannot be given in the general case, checking this explicitly at 
each extremum is also problematic. One can circumvent this computational difficulty using 
the following Proposition: 
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Proposition 1. For each extremum, the hessian remains signature- definite within each conical 
region described in the previous subsection. 

Proof. Let us fix and move continuously in the parameter space, keeping it strictly 
inside one of the regions described in the previous subsection. Let us pick up an extremum and 
follow how it changes when A 11 moves. Its position, its depth as well as the eigenfrequencies 
are algebraic functions of the components of A^ and therefore they also change continuously. 
So, if the hessian changes signature along at the endpoints of some A^ path, then there exists 
a point, at which one of the eigenfrequencies is zero. Thus, the expansion of the potential at 
this point starts from the third or fourth order term, and this points corresponds to merging 
of two or three simple extrema, respectively. 

But such a merging cannot happen for A^ lying strictly inside the caustic regions. Indeed, 
if two sufficiently close points r^ 1 and are both extrema of the potential, then their respective 
zeroth components r a0 and r b0 are also close, so, the intersection points of the corresponding 
ellipsoids lie close to the boundary of a caustic region. In the limit r% — > r^, the intersection 
points, An being among them, approach the envelope (loosely speaking, the envelope can be 
viewed as the locus of intersections of the "successive" ellipsoids). □ 

So, since the signature of the hessian remains the same for all A^ inside some region, 
one can select some representative A^, calculate the signature of the hessian for it, and then 
extrapolate the results for all the points inside this region. 

Let us now calculate the number of minima inside the innermost region of the A^ space, 
see Fig. (TJ right. For the representative point in this region, A^ = (A , 0, 0, 0), calculations 
can be easily done explicitly. Indeed, it follows from ( fl8l) that there are three pairs of extrema 
at r = Tq , , and 7q : 

rS 1} (l, ±1,0,0), rf(l, 0, ±1, 0), rf(l, 0, 0, ±1) . (23) 

Again, let us order the eigenvalues B i: B\ < B 2 < B 3 and expand the potential near the point 
(r M } + = r^\l, 0, 0, 1). If r M = r (l, sin^cos^, sin 8 sin 0, cos6) = (r M ) + + 5r M , then 

V = -A r + K 2 (B - B x sin 2 6 cos 2 <p - B 2 sin 2 9 sin 2 <p - B 3 cos 2 6) 

Al | B -B 3 ^ )2 | A 2 g 2 (j? 3 -i?i)cos 2 0+(i?3-i? 2 )sin 2 ^ 



2(^0-^3) 2 v uy 2 OBo-Bg 



Here, 5ro and 9 are small, while can be arbitrary. Since B 3 is the largest spacelike eigenvalue, 
this point is a minimum, and so is the other extremum of this pair, (r")_ = rj 3) (l, 0, 0, -1). 
The same calculation for the extrema at r = and r^ 2 \ shows that they are saddle points. 
Thus, we find that for A^ lying in the innermost region, the potential has two separate minima 
and four separate saddle points in the orbit space. 

As A^ moves out of this region, the number of minima does not increase. Indeed, one 
can show that crossing the caustic surface at a generic point leads to disappearance of two 
saddle points or of one saddle point and one minimum, but it cannot, for example, lead to 
disappearance of three saddle points and appearance of a new minimum. This can be also 
verified with the straightforward calculations similar to (12^1) by selecting points A^ lying on 
the axes (for this choice, all the extrema can be also studied with explicit algebra). Therefore, 
we arrive at the following Proposition: 



12 



Proposition 2. The most general quadratic plus quartic potential with two order parameters 
can have at most two distinct local minima in the orbit space. 



3.3 The principal caustic cone 




Figure 2: The same as in Fig. [TJ right, but with number of minima and saddle points shown 
separately as N minima + N sadd[e . Thick lines show the section of the principal caustic cone. 

In Section 13.11 we showed that for a generic potential there exist two caustic cones in the 

it) (2) (2) (3) 

space. If Bi < B 2 < B3, they corresponds to r < r < r and r < ro < r^ , respectively. 
The analysis of Section 13.21 shows that these two caustic cones play different role. It is the 

(2) (3) 

second cone, with Tq < ro < r , which we call the principal caustic cone, that separates 
regions with different number of minima. This is illustrated by Fig. [2] which is an updated 
version of Fig. [U right, with the numbers of minima and saddle points shown separately. The 
most straightforward proof is based on the stability analysis of the extrema in three situations 
with lying on each of the three axes and Proposition [IJ The other cone just separates 
regions with different numbers of saddle points and does not affect directly the properties of 
the global minimum. 

3.4 Geometric reformulation of the search for the global minimum 

Consider again the potential term in (fi"3l : 

V = -A^ + ~B lxu r»r v . (25) 

Let us exploit the freedom in definition of B^ u to make B > and all Bi < 0. Then (B^ 1 )^ 
exists, and fT23|) can be rewritten as 

V^B^-aPyf-a^ + Vo, a^iB- 1 )^, V = ~{B- 1 )^A V . (26) 

Let us now define an equipotential surface Ai c as a set of all vectors p^ in the Minkowski 
space M. such that 

B^f = B pl + \Bi\vl = C ■ (27) 
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One sees that equipotential surfaces exist for C > and are 3-ellipsoids nested into each other, 
with their eigenaxes aligned in the .B^-diagonal frame along the eigenaxes of B^. 

Returning to the potential fl26|) we see that C is related to the values of the potential: 
C = 2{V — Vq). Therefore finding points in the orbit space with the same value of V amounts to 
finding intersections of the corresponding Ai° with the forward lightcone LC + . In particular, 
to find a local minimum of the potential in the orbit space, one has to find an equipotential 
surface that touches LC + (we say that two surfaces "touch" if they have parallel normals at 
the intersection points). The glo bal minimum corresponds to the unique equipotential surface 
A4 Cmin that only touches but never intersects LC + . 




Figure 3: A 1 + 2-dimensional illustration of the contact between A4 mln and LC + . Shown is 
the case of a M lying outside LC + . 

Thus, the geometric strategy for the minimization of the potential is the following: 

• Construct a family of 3-ellipsoids Ai c at the base point a M ; 

• Find the unique 3-ellipsoid A4 min that merely touches LC + but never intersects it; 

• The contact point or points give the values of r M (hence, ipi and 1P2) that minimize the 
potential. 

To facilitate the visualization, Fig.[3]shows a 1+ 2-dimensional analogue of the contact between 
Ai Cinin and LC + . In this particular example, a M , which is located at the center of the ellipsoid, 
lies outside LC + . 

Alternatively, using (|27|) one can interpret the potential as the distance squared from the 
point d 1 in the Euclidean metric diag(i?o, I-B2I, \Bs\). The minimization problem is then 
reformulated as a search for points on LC + that are closest to a M in this metric. Since the 
forward lightcone LC + (together with its interior) is a concave region, this representation 
immediate leads to the following conclusion: the necessary condition for the existence of a 
degenerate minimum is that a M lies inside LC + : a^a^ > 0. Later, in Section [5731 we will give 
necessary and sufficient conditions for this to happen. 

The two geometric constructions described above, the ones based on the equipotential 
surfaces and on the caustic cones, are related to each other in the same manner as a planar 
curve to its evolute. To illustrate this relation, let us consider a simple planar problem: find 
on the unit circle points of local minima of the "potential" V = Bij(rii — ai){rij — aj) with 
Bij = diag(5i, B 2 ), B\ ^ B 2 , Bi > 0. In the coordinates fii = (\^B)ij7ij, the "potential" 
becomes V = \hi — ai\ 2 , while the unit circle is transformed into an ellipse. One can easily 
verify that the number of local minima depends on whether point a 1 lies inside the evolute of 
this ellipse. 
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4 Potentials stable in a weak sense 



Let us now discuss how the above constructions change for a potential stable in a weak sense, 
i.e. a potential whose V4 can have flat directions in the ^j-space, along which the potential is 
stabilized by the V 2 term. 

A flat direction of the quartic part of the potential in the ^-space corresponds to a vector 
r M G LC + in the orbit space such that B^r^r" = 0. Such an r M must be an eigenvector of B^, 
see Appendix [B] Let us first assume that there is only one such direction. Aligning it with 
the first axis, one can diagonalize B^ v in the "transverse space" , bringing it to the following 
generic form: 



+ SB 


SB 







5B 


-B + 5B 










-B 2 













-B 3 


) 



with SB > 0, B > B 2 , B*. (28 



( B ( 
V 

Note that in contrast to the potentials stable in a strong sense, ( ITBI) . this B^ v cannot be 
diagonalized by an 50(1, 3) transformation. Indeed, a boost along the first axis with "rapidity" 
rj leads to the same B^ v as ( 1281) but with redefined SB — > e~ 2r) SB (see Appendix [A} . If SB ^ 0, 
then B^ v is not diagonalizable. 

To find the number of extrema in this case, one can again start with system (fl6l) but instead 
of considering fixed Aq sections in the A^ space one can fix one of its lightcone components. 
Let us introduce the lightcone decomposition of any four-vector: 

p"=p + 7i£+p_T£+p£, < = (1, ±1, 0, 0), (29) 

where = (0, 0, p 2 , P3). The lightcone coordinates p± are related to the zeroth and first 
coordinates po, pi as p± = (po ±pi)/2. Then, system (fl6l) can be rewritten as 

(B - A)r_ = A_ , 
(B -X)r + + 2SBr^=A + , 
(B i -\)r i = A i , 2 = 2,3, 
Ar + r^ = r\ + r| . 

The condition that the quadratic part of the potential V 2 = —A^r^ increases along the flat 
direction of V4 given by n+ implies that A^ < 0. Repeating the geometric analysis described 
in full detail for the potential stable in a strong sense, we arrive at the following conclusion: 



If A^ lies inside the past lightcone LC , then system ( 1301) has no solution. The global 
minimum is at the origin. 

If A^ lies outside LC~ (but still with A_ < 0), then system (13T)j) has a unique solution. 
The corresponding unique extremum of the potential is the global minimum. 



Suppose now that there are more than one flat direction of V4. Let us pick up two such 

B^ v r\ = Ai7"i M , B^ v r v 2 = X 2 r 2f , . 



distinct vectors rf, G LC + , both eigenvectors of B 
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Since r^r 2fl 7^ 0, one obtains Ai = A 2 . Then, using the freedom in definition of B^ u , one 
can always set this common eigenvalue to zero. Then for any linear combination of these two 
vectors one gets 

B P (on + (3r 2 ) v = . (31) 

Consider now such an 5*0(1, 3) transformation that makes oc n+ and r% oc nt. Then, B^ v 
that satisfies (13ip takes the following generic form: 



B 



( 

































-B 2 







Vo 








-B 3 


J 



with > B 2 , B 3 



(32) 



This form is diagonal (note also that it is equivalent to (128]) with 5B = 0), so one can again 
switch to the fixed A sections. Since — A M r M must increase along n±, A < 0. Repeating the 
same analysis as in Section [3TT| one finds that the potential has a unique non-trivial extremum 
if I Ai| < \A Q \ and A^A^ < 0. 

So, a potential stable in a weak sense is similar to the potentials stable in a strong sense 
with Aq < in the frame with diagonal B^ u . It can have no more than one non-trivial 
extremum, which is then necessarily the global minimum. 



5 Symmetries and their violation 
5.1 Explicit symmetries 

As explained in Section 12.11 the free energy remains invariant under an appropriate simulta- 
neous transformation of the order parameters ipi and the coefficients. It can happen, however, 
that the free energy is invariant under some specific transformation of ipi (or the coefficients) 
alone. We call this symmetry an explicit symmetry of the free energy. 

In the orbit space, this symmetry corresponds to such a map of the Minkowski space M. 
that leaves invariant, separately, B^ v r^r v , A^r^ and K^p^. The notion of explicit symmetry 
is invariant under the Lorentz group of the orbit space transformations; so, any SO(l,3) 
transformation leaves a given model in the same symmetry class. 

In simple cases the presence of a symmetry can be evident from a direct inspection of the 
free energy functional, see e.g. examples in Section [71 In fact, in many concrete applications, 
the Landau potential is even constructed in such a way that some symmetry is preserved. In 
the general case, however, a non-evident hidden symmetry can exist even in complicated forms 
of the free energy, without being easily noticeable. So, one needs a reparametrization-invariant 
criterion that can help recognize the presence of a symmetry using only K^, A^, and B^ v . In 
addition, it would also useful to know what this symmetry is. 

Both questions are answered by the following Proposition: 

Proposition 3. Suppose that the free energy ^ is explicitly invariant under some transfor- 
mations ofr^. Let G be the maximal group of such transformations. Then: 

(a) G is non-trivial if and only if there exists an eigenvector of B^ u orthogonal both to A^ and 

i< lt ; 

(b) group G is one of the following groups: Z 2 , (Z 2 ) 2 , (Z 2 ) 3 , 0(2), 0(2) x Z 2 , or 0(3), and 
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depends on the number of the eigenvectors of to which and are orthogonal, and on 
whether B^ u has degenerate eigenvalues. 

Proof. Let us start with the potential stable in a strong sense. Let us denote the group of all 
explicit symmetries of B^ v , A^, and by Gb, Ga, and Gk, respectively. Obviously, 

G = G B n G A fl G K . (33) 

The group of explicit symmetries is necessarily a subgroup of the 0(3) transformation group 
of the 3-dimensional space in the -B Mi ,-diagonal frame; so one can switch to the spacelike parts 
only (Bij, A h Ki). 

Consider now Gb- If all spacelike eigenvalues of B^ are different, then its only symmetries 
are reflections of each of the spacelike eigenaxes, which generate group Gb = {Z2Y ■ If two 
eigenvalues coincide, then Gb = 0(2) x Z 2 , and if all three of them are equal, then Gb = 0(3). 
Note that in all of these cases the following statement holds: if some Z 2 group is a subgroup 
of G^, then its generator flips the direction of some eigenvector of B^. 

Similarly, Ga is 0(2) (rotations around the axis defined by Ai), if Ai is a non-zero vector, 
and 0(3) otherwise. The same holds for Ki, the only difference being the direction of the 
axis. If we want G to be non-trivial, then the lowest possible symmetry of Ai and Ki together 
(given by a Z 2 group) must be also a symmetry of Bij, i.e. it must flip one of the eigenvectors 
of B^. In other words, both Ai and Ki are orthogonal to this eigenvector. Since this purely 
spacelike eigenvector is also the eigenvector of B^ u , we arrive at the first statement of this 
Proposition. 

Detailed classification depends on the number of eigenvectors of B^ that are orthogonal to 
Ai and Ki. 

• If Ai and Ki are orthogonal to all three eigenvectors (Ai = Ki = 0), then G = Gb- 

• If Ai and Ki are orthogonal to two eigenvectors (y4j||ii"j and are themselves eigenvectors 
of B^), then G = (Z 2 f or 0(2). 

• Finally, if there is only one eigenvector of B^ orthogonal both to Ai and Ki, then the 
symmetry group is Z 2 . 

For a potential stable in a weak sense, we first note that the eigenvectors of B^ v are either 
the lightcone vectors or purely spacelike eigenvectors. Since K^ lies inside LC + , in cannot be 
orthogonal to any lightcone vector. Therefore, one has to check the above conditions only for 
the "transverse" eigenvectors, which reduces the above list of possible symmetry groups to Z 2 , 
(Z 2 f,0(2). ' □ 

The necessary and sufficient condition formulated in the first part of this Proposition can 
be written in a reparametrization-invariant way. The method is essentially the same as in 
[3TI [32] and is based on a simple observation: if a 3- vector Oj is orthogonal to some eigenvector 
of a real symmetric matrix by, then the triple scalar product of vectors a«, hjOj, and hjbjkak 
is zero. In Minkowski space we introduce 

K^ = K^, K X p = B»K v , K 2 ^ = (B 2 )^K U , K 3fl = (£? 3 )/ K v , (34) 
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where B k is the k-th power of B^ v . The same series can be written for A^. For any four 
four- vectors a' 4 , 6 M , c M , and we introduce the short-hand notation 



(a, b, c, d) = e^ pa a^b u c p d a . 

Then the condition "there exists an eigenvector of B pv orthogonal to K p can be written as 

(K ,K h K 2} K 3 ) = 0. (35) 

Note that since K p always lies inside the future lightcone, it can be orthogonal only to spacelike 
eigenvectors of B^ u , which is exactly what is needed. Then, the statement of Proposition [3^ 
can be reproduced if we accompany fl35l) with a similar condition for A^, 

(Ao,A 1 ,A 2 ,A 3 )=0, (36) 

and with the condition that these two 4- vectors be orthogonal to the same eigenvector of B pvi 
for example: 

(A ,A 1 ,A 2 ,K ) = 0, (37) 

Conditions (j3"5l) - (l3~Tj) can be straightforwardly checked in any basis once B^ u , A^, and K M 
are known. Thus, the presence of any hidden symmetry can be verified without the need to 
find this symmetry explicitly. 



5.2 Symmetries of the potential vs. symmetries of the free energy 

Explicit symmetries of the entire free energy depend on By, Ai, and Ki, while the symmetries 
of the potential depend only on Bij and A^. Therefore, it might happen that the potential has 
a larger symmetry group than the entire free energy. A simple example is 

F = k (\D^\ 2 + |^ 2 | 2 ) + 16A (|^| 2 - v 2 ) 2 + A (> 2 | 2 - Av 2 ) 2 . (38) 

The potential here is symmetric under ip2 *-* 2^i, while the gradient term is not. 

The two notions, i.e. the symmetry of the potential or of the entire free energy, play 
different roles. When one seeks for the minimum of the Landau potential, the coefficients 
in the gradient term (K^) are irrelevant. However, the symmetry of the spectrum of small 
oscillations of the order parameters above the ground state is the one of the entire free energy 
functional. 



5.3 Spontaneous breaking of an explicit symmetry 

Even if the Landau potential is invariant under some transformation of $, the values (<&) that 
minimize it do not necessarily have to preserve the same symmetry. In the orbit space, if the 
Landau potential is invariant under a group G of transformation of r M , then the position of 
the global minimum might be invariant only under a proper subgroup of G. In such situations 
one speaks of spontaneous breaking of the symmetry. Since the set of all global minima is 
invariant under the full explicit symmetry group G, the spontaneous breaking of an explicit 
symmetry always leads to degenerate global minima. 

For our problem, several results follow immediately from Proposition [2J 
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1. The global minimum can be only twice degenerate. 



2. Minima that preserve and violate any discrete symmetry cannot coexist. 



3. The maximal breaking of any discrete symmetry consists in removing only one Z 2 factor: 
7 nf^(Z 2 ) 



:z 2 ) k ^ (Zo) k -\ with A; = 1,2, 3. 



In addition, in [29] it was proved that the twice degenerate global minimum of Landau 
potential with quadratic and quartic terms is always realized via spontaneous breaking of some 
explicit Z 2 symmetry of the potential (but not necessarily of the entire free energy!). 

Let us now consider the question when a given explicit symmetry is broken, focusing on 
the discrete symmetry case. 

First of all, the global minimum must be degenerate. This immediately leads to the 
conclusion that the spontaneous violation can take place only in potentials stable in a strong 
sense, and in addition, only when lies inside the principal caustic cone. To make the 
discussion concrete, consider A^ and in the 5^-diagonal frame. Suppose that all Bi are 
distinct and the components A 3 = K 3 = 0, while the other components are non-zero. Then, 
the free energy has an explicit Z 2 symmetry generated by reflections of the third coordinate. 
This explicit symmetry is conserved, if the global minimum is at r M = (r , r±, r 2 , 0), and it is 
spontaneously broken if the two degenerate global minima are at r± = (ro, n, r 2 , ±7*3) with 
r 3 ^ 0. 

Let us now recall the "shrinking ellipsoid" construction of Section 13.11 A degenerate 
extremum implies that two distinct points Ui±, when inserted in system ffl8|) . give the same 
point Ai = (Ax, A 2 ,0) at the same r . This happens only when r = Tq and the planar 
vector (Ax, A 2 ) lies inside the ellipse with semiaxes 



\Bx — B 3 \ \B 2 — B-, 



A ^ ^, A - 



3 1 



Bq — B 3 Bq — B 3 

Besides, as we prove in Appendix [Cj in order for this extremum to be minimum, B 3 must be 
the largest (i.e. the closest to Bq) eigenvalue among all Bi. Thus, one arrives at the following 
necessary and sufficient reparametrization-invariant criterion for the spontaneous violation of 
a Z 2 symmetry (along the third axis): 



A 2 A 2 A 2 



Bs>Bx,B 2 and - 2 D , 2 < /D " , 2 ■ (39) 



(B 3 -Bx) 2 (B 3 -B 2 ) 2 (B -B : 
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It immediately follows from here that a M defined in Section [3^1 lies inside the forward lightcone. 

Finally, we would like to stress one important point. What is relevant for the whole 
discussion is the group G of explicit symmetries and its reduction upon symmetry breaking, but 
not the particular realization of the transformations of this group. For example, the free energy 
can be symmetric under ipi — > ip* transformation or under ipx ^2 transformation. These 
seemingly distinct Z 2 symmetries are, in fact, just different realizations of the same symmetry 
class. This becomes evident in the orbit space, as the former transformation corresponds to 
the flip of the second axis, while the latter one corresponds to the flip of the third axis; so, 
both models are related by a reparametrization transformation. Therefore, all properties of 
spontaneous violation of these two particular sorts of the Z 2 symmetry are in fact identical. 
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Thus, there exists a reparametrization-invariant class of ^-symmetric models, a reparametrization- 
invariant class of (Z 2 ) 2 -symmetric models, etc. Our analysis applies to all particular realiza- 
tions of a given symmetry in contrast. This rather simple fact illustrates the usefulness of 
considering the most general GL model and might lead to establishing of direct links between 
seemingly unrelated models. 

6 Phase diagram and phase transitions 
6.1 Phase diagram 

The results obtained in the previous sections allow us to explicitly describe the phase diagram 
of the most general model with two order parameters in the mean-field approximation. We do 
this by classifying the phases according to the symmetries of the free energy functional and 
to the properties of the ground state. For definiteness, we again sort the eigenvalues Bi as 
B\ < B 2 < B 3 . 

• Potential is stabilized by the quadratic term (A 11 lies inside LC~). The global minimum 
is at the origin, ipi = ip 2 — 0; this is the high-symmetry phase. 

• Potential stable in a weak sense. The global minimum is always non-degenerate and 
preserves any explicit symmetry of the free energy. 

• Potential stable in a strong sense. B^ v can be diagonalized, and one can work with its 
spacelike part. 

— All Bi are distinct. 

* Ai and Ki are generic vectors (not orthogonal to any eigenvector of By). No 
explicit symmetry is present. There can be one or two non-degenerate minima, 
depending whether Ai lies inside the principal caustic region. 

* Ai and Ki are both orthogonal to one eigenvector of Bij. The explicit symmetry 
group is Z 2 . The ground state can either break or preserve this symmetry. The 
symmetry is broken if it is the third axis that Ai and Ki are orthogonal to (i.e. 
A 3 = K 3 = 0) and if condition fl39|) is satisfied. 

* Ai and Ki are both parallel to the same eigenvector of B^. The explicit sym- 
metry group is (Z 2 ) 2 . The ground state can either preserve this symmetry or 
break it to Z 2 . The criterion of the symmetry breaking is the same, (1391) . apart 
from the fact that now one of A\, A 2 is zero. 

* Ai and Ki are orthogonal to all three eigenvectors of B^ (i.e. Ai — Ki — 0, 
A > 0). The explicit symmetry group is (Z 2 ) 3 . The global minimum is always 
twice degenerate and breaks this symmetry to (Z 2 ) 2 . 

— Two eigenvalues among Bi coincide. Case B\ < B 2 = B 3 . 

The principal caustic cone reduces to a segment along the first axis. 

* Ai is not aligned along the first axis. Then, the global minimum is non- 
degenerate. If Ki lies in the (emj, Ai) plane, then there is an explicit Z 2 
symmetry which is preserved at the minimum. 
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* Ai is the eigenvector of along the first axis (A 2 = A 3 = 0). Then, if 
condition fl39|) is satisfied, there is a continuum (namely, a circle) of degenerate 
minima, otherwise, the global minimum is non-degenerate. If, in addition 
K 2 = K 3 = 0, then the explicit symmetry group is 0(2), and in the case 
of symmetry breaking, it is broken to Z 2 . Instead, if Ki is a generic vector, 
then the explicit symmetry group is Z 2l and it may be preserved or broken 
depending on which minimum in the continuum is selected. 

* Ai = Ki = 0, Aq > 0. The explicit symmetry group is 0(2) x Z 2 . There is 
always a continuum (a circle) of global minima, and the symmetry is broken 
to(Z 2 ) 2 . 

— Two eigenvalues among Bi coincide. Case B\ = B 2 < B 3 . 

The analysis is similar to the case with all distinct B i7 with the following differences: 

* If Ai and Ki are both parallel to the eigenvector e^, the explicit symmetry 
group is 0(2) and it is always preserved at the global minimum. 

* If Ai = Ki = 0, then the explicit symmetry group is 0(2) x Z 2 , which is broken 
at the global minimum to 0(2). 

— All three eigenvalues Bi coincide. 

* For a generic pair of Ai and K iy there is always an explicit Z 2 symmetry. It 
is promoted to the 0(2) symmetry if Aj and Ki are parallel, and to the 0(3) 
symmetry if Ai = Ki = 0. The symmetry is always preserved, apart form the 
case Aq = 0, when it is broken to 0(2) or Z 2 . 



6.2 First and second order phase transitions 

A remarkable property of the two-order-parameter model is that it can have a first-order 
phase transition even at zero temperature and in the mean-field approximation. It is entirely 
due to coexistence of two local minima (in the orbit space) with different depths. If upon 
continuous change of the coefficients the relative depth of the two distinct minima changes sign 
(the shallower minimum becoming the deeper one), the system occupying initially the global 
minimum becomes metastable and can jump into the new global minimum via fluctuations or 
quantum tunneling. 

If B 3 is non-degenerate, then the surface of first-order phase transitions in the Ai space is 
the interior of the ellipse at i.e. it is given by Ai = (Ax, A 2 , 0), where Ax, A 2 satisfy (|39|) . 
The border of this ellipse, 

A 2 A 2 A 2 



A 3 = , — i— - + 



' (B 3 -Bx) 2 (B 3 -B 2 ) 2 (B -B 3 ) 2 ' 

is the critical line, at which second-order phase transition takes place. If B 3 is degenerate, 
B\ < B 2 = B 3 , then the second order phase transition takes place at isolated points 

B 3 — Bi 



A 3 = A 2 = , Ax = ±\ 



Bo — B 3 ' 



and the points of the first-order transitions form a linear segment between them. Finally, if 
all Bi are degenerate, then there is a single critical point at the origin, Ai = 0, and there is 
no first order phase transition. 

Reconstruction of critical surfaces/lines in the A^-sp&ce is obvious. 
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6.3 Critical properties: an example 

It appears plausible that all the mean-field critical exponents of the general 20P GL model 
are of geometric nature and can be calculated without the knowledge of the exact position of 
the global minimum. Here, we do not explore this issue in full detail, but just illustrate it 
with one example: calculation of the critical exponent of the correlation length with the aid 
of differential geometry. 



generic critical 



critical 




Figure 4: (a) Schematic view of the principal caustic region. The two arrows enter it via a 
generic or a critical point, (b) Schematic change of the potential upon a generic entrance into 
the principal caustic cone, (c) Schematic change upon the entrance via a critical point. 

Let us fix B^ u and move continuously in its parameter space. When it crosses the 
principal caustic cone (see Section 13.31) . the number of minima changes. Two scenarios are 
possible, Fig. Hh- If A^ 1 enters the principal caustic cone through a generic point, then an 
additional local minimum appears together with an additional saddle point, as it is shown 
schematically in Fig. [4b. This bifurcation does not involve the global minimum. However, if 
A^ enters the principal caustic cone through any of the critical points, then it is the global 
minimum that bifurcates into minimum-saddle-minimum sequence, see Fig. [4b. 

When A^ approaches the critical surface/line, the corresponding eigenfrequency decreases 
and turns zero exactly at the critical surface. If the distance from A4 to the critical surface is 
e — > 0, the eigenfrequency associated with this bifurcation decreases as 

uj 2 oc e 5 . (40) 

Here we used the fact that the Jacobian of the map of non-Goldstone modes in -i/^-space to the 
surface of LC + is regular, if the extremum is not at the origin, see Eq. (122]) . The correlation 
length then behaves as r c oc e~ s / 2 . We argue that the value of S is of geometric nature and 
can be calculated without the knowledge of the exact position of the minimum. 

Let us first note that in the case of a single order parameter 5=1 simply because the 
eigenfrequency uj 2 is linearly proportional to the coefficient a in the potential ([1]). In the two- 
order-parameter case, due to the higher dimensionality of the A M -space, one can approach a 
critical point from different directions. 

As we described in Section 13.41 the search for the global minimum can be reformulated 
as a search for points lying on LC + that are closest to a given point a M in the Euclidean 
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1 



Figure 5: (Color online) The principal caustic region in the Aj space (dashed line) and its 
border, the critical line (think solid line). The critical exponent 5 depends on whether 
approaches the critical line always staying in the plane A 3 = or from outside the plane. 

metric diag(i?o, [-Si J , I-B2I) | -B3 1 ) • In Appendix IDl using a planar example, we show how to 
apply differential geometry to analyze the properties of the potential near a critical point. We 
showed, in particular, that the value of 5 depends on the direction of approach to the critical 
point: 

generic direction — > 5 = 2/3, 
symmetric approach — > 5 = 1, (41) 

These exponents are robust in the sense that they remain the same for almost all, in the 
measure-theoretic meaning, regular planar curves. It applies also to the second order curves 
which share the key property of the generic curves that they have points of no more than 
4th-order contact with a circle. 

This technique can be extended to higher dimensions leading to the same results. So, 
exponents fj4T|) apply to our problem as well, where "symmetric approach" is understood as 
11 A4 lying in the A 3 = plane" , see Fig. [5j 

It would be interesting to check the critical properties of all possible phase transition in 
20P GL model and see how many classes of critical behavior it can incorporate. 

7 Examples 

Here we illustrate the general approach with several simple examples. Some of them are 
relevant for condensed- matter problems discussed in literature. 

7.1 Real coefficients 

Consider a free energy functional ([!]) with all real coefficients. It implies that the free energy 
remains invariant under simultaneous transformation if)^ — > if)*, which corresponds in the orbit 
space to the reflection of the second coordinate: 

r » = (r , n, r 2 , r 3 ) -> (r , n, -r 2 , r 3 ) . (42) 
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Consequently, the four-vectors and four-tensor of the coefficients are: 



B 



hi/ 







V 






-B 2 




o, ■) 



(43) 



/ 



where dots indicate generic values. Evidently, conditions ( l35l) - (|371) are satisfied. 

Just to give a particular example, consider the free energy functional of a two-gap su- 
perconductor in the dirty limit, see Eq. (49) in [14] . Its potential can be rewritten in the 
reparametrization-invariant way ffl~3]) with 



B 



-2bi 


bi-&2 
















bi- 






b 1 +b 2 



A" = - 
2 



v 2 2 

cli — a 2 , di, 0, ai — a 2 ) . 



h ) 



(44) 



Here we used the notation of [14] : coefficients a±, b\ and a 2 , & 2 refer to the properties of the 
first and second order parameters, respectively, while a^, hi describe interaction terms. The 
gradient terms considered in [H] are anisotropic, but they also contain real coefficients. As a 
result, A 11 and K^ 1 are orthogonal to the second eigenvector of B^ u . Note that in this example, 
the eigenvalue B 2 = 0. In order to find the other eigenvalues, one has to solve the cubic 
characteristic equation. 

If the position of the global minimum has (r 2 ) = 0, then the symmetry is preserved, and 
there is no relative phase between the two order parameters. If (r 2 ) ^ 0, then the symmetry is 
spontaneously broken, and (tpi) and (?/> 2 ) in the ground state have a non-zero relative phase. 
In order to find whether spontaneous violation takes place, one has to diagonalize B^, find 
its eigenvalues as well as find A^ 1 in this frame, and then check inequality (1391 . 




Figure 6: Phases of the general Ginzburg-Landau model with two real order parameters on the 
Ai plane, classified according to the number of minima and conservation/violation of explicit 
symmetries. Shown is the case Bi < B 3 . 
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If one requires, in addition, that the order parameters ipi themselves be real, then r 2 = 0, 
and the second axis can be omitted altogether. The orbit space is then simplified to the forward 
lightcone in the 1 + 2-dimensional Minkowski space. Repeating the analysis of Section [3], one 
obtains now at most four non-trivial extrema, among which up to two can be minima, and 
only one caustic cone. 

The phase diagram in this case is simpler. For example, in Fig. [6] we show phases in the 
Ai space for the case B 3 > B\ and, for simplicity, we assume that Ki — 0. The astroid shown 
here is the planar analogue of the cusped region from Fig. [I] and Fig. [2j It separates the 
Ai regions corresponding to potentials with one minimum (phase I) and two minima (phase 
II). In addition, if Ai lies on the axes, the potental has an explicit symmetry. Dashed lines 
corresponds to the cases when the ground state conserves the symmetry (phases Ic and lie), 
while the thick solid line corresponds to the phase that spontaneously break the discrete 
symmetry (phase IIv). 

7.2 No ifji/i/j2 mixing in the quartic potential: 65 = 6g = 

The situation simplifies considerably if the quartic potential does not mix ipi and ifi 2 , i-e. when 
^5 — ^6 = 0. In this case B^ u breaks into two 2x2 blocks and can be easily diagonalized by 
a boost along the third axis (see Appendix |X] for details) with "rapidity" 

where we assumed b\ > b 2 , and by the rotation between the first and second axes by an angle 
equal to the half of the phase of 64. The resulting eigenvalues are: 

B = \{\fhb 2 + 63) , B 1>2 = ±^, B 3 = ^{-y^b~ 2 + b 3 ). (46) 
The condition for the stability in a strong sense of the potential is 

61 >0, b 2 >0, Vlhb~ 2 + b 3 > \b 4 \ . (47) 
If A^ and are generic vectors, then the further analysis proceeds as in the generic case. 

7.3 Interaction only via the l^il 2 ^! 2 term 

Let us assume now that the only interaction between the two OPs is given by the £>3 1 1 2 1 "V^s 1 2 
term. The four- vector A^ can be written as (A , 0, 0, A 3 ), so that the potential has an explicit 
0(2) symmetry. The eigenvalues of B^ v expressed in terms of the original coefficients are 

B = \{Vhh + 63) , £1,2 = 0, J B 3 = i(-v / M 2 " + &3), (48) 

where for stability we require \/bib 2 + 63 > 0. In the frame where B^ v is diagonal, A^ takes 
the form 

(A , 0, 0, A 3 ) = (A cosh i] + A 3 sinh rj, 0, 0, A 3 cosh 77 + A sinh rj) , (49) 
where 7] is given by ( )45l) . 
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If B 3 > 0, i.e. 63 > \/M>2; then the explicit 0(2) symmetry is always conserved, since Bx 
and B 2 are not the largest spacelike eigenvalues. The global minimum is at (r M ) oc n+ or nt, 
which corresponds to 



(V> 1 )^o,(v 2 ) = o. 



or 



<Vi) = 0,(^)^0. 



If 53 < and if A' 4 lies inside the caustic cone, then this symmetry is spontaneously broken, 
and there exists a continuum of global minima with both (ipx) 7^ and (^2) 7^ and an 
arbitrary relative phase between them. The condition that A 11 lies inside the caustic cone is 



A, 



B\ — B 3 



< 



A, 



B — B 1 



which, in terms of the original coefficients translates into 

ax— < a 2 < a\ — . 
bx h 



(50) 



Of course, the same bounds can be obtained by direct calculations. 

Note that in the case of no interaction at all, 63 = 0, the condition for the symmetry 
violation reads \A 3 \ < A . It means a\ > 0, a 2 > 0, which is indeed expected. 



8 Solitons 

Two local order parameters can lead to existence of solitons, i.e. states with non-trivial co- 
ordinate dependence of the mean- field values of the order parameters (0j)(r) stable against 
small variations 5{ipi)(r) of these OP profiles. Some particular versions of such solitons have 
been already described in literature. For example, in [33], a one- dimensional two-band su- 
perconductor with a simple interband interaction term was considered, whose ground state 
corresponded to (^1) and (^2) with zero relative phase. Then, a typical sine-Gordon soliton 
was constructed with the relative phase between the two OPs continuously changing from 
zero to 2tt at x = ±00, correspondingly. Similar solitons in the scalar sector of 2HDM were 
described in [M]. 

Existence of solitons in a given 20P GL model depends on the geometry of the potential in 
the orbit space. For example, in order to support a one-dimensional soliton similar to the one 
described above, the Landau potential must have a certain "valley" (i.e. a region of low values 
of the potential) of non-trivial topology on the forward lightcone LC + , that would include the 
global minimum and a saddle point. At x — > —00, (ipj)(x) approach their values at the global 
minimum. As x increases, the corresponding point (r M ) (x) in the orbit space moves away from 
the global minimum position, follows some path in the valley and returns again to the global 
minimum. Small variations of (ipi)(x) would pull this path out of the valley, increasing its 
potential energy. 

Existence, stability and geometric properties (e.g. dimensionality) of these solitons are 
sensitive only to the general structure of the model, and do not require one to search for the 
explicit position of the extrema. Therefore, one can hope to obtain these criteria for a general 
20P GL model in terms of geometric constructions studied in this paper. 
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9 Multicomponent order parameters 



So far, we assumed that the order parameters ipi and ^2 are just complex numbers. However, in 
many physical situations one introduces multicomponent order parameters. Examples include 
2HDM, superfluidity in 3 He, non-conventional superconductivity, spin-density waves, etc. 

The formalism presented above is applicable to these cases as well. In fact, it was first 
developed in [2U|29j specifically for 2HDM. Here, we discuss characteristic features that appear 
in a generic GL model with two iV-vector order parameters and a [/(iV)-symmetric potential. 



9.1 Modifications to the formalism 

Let us assume that each ipi is an iV-dimensional complex vector: ip ia , a = I, . . . , N. A 
L r (A^)-symmetric potential must depend on the order parameters only via scalar combinations 
(ipjipj), i,j,= 1,2, which parametrize the [/(iV)-orbits. The only difference with the scalar 
case is that an additional term proportional to 

{4M(4^i) ± l^ilV 2 | 2 



(51) 



appears in the potential, with a new independent coefficient 63 in front. The definition of r M 
remains the same; however 



> 0. 



(52) 



Therefore, the orbit space now is not only the surface, but also the interior of the forward 
lightcone LC + . This removes the degree freedom in definition of B^, making it uniquely 
defined: 
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(53) 



The requirement that the potential is stable in a strong sense implies that B^ v must be positive 
definite on and inside LC + . This leads not only to B > B iy but also to B > 0. Note that 
due to the absence of freedom in B^ u , cases with singular B^ v and with Bi of different signs 
must be considered as well. 



9.2 Consequences 

Let us discuss the modification of the above analysis due to the multicomponent order param- 
eters. 

A new phase appears, which is characterized by a stronger breaking of the symmetry of the 
potential. It corresponds to the global minimum of the potential (r M ) lying strictly inside the 
lightcone LC + . This is possible only when and (^2) are not proportional to each other. 
In other words, one can always perform a simultaneous "intra- vector" U(N) transformation 
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of both order parameters that makes them 



( o \ 



V«i J 



( o \ 



<^2> 



u 



(54) 



where dots indicate zeros. Here, u, vi, t>2, £ are real, and u and Vi must be non-zero in order 
for r M constructed from them to lie strictly inside LC + . 

Solution (}54"j) with non-zero u preserves only a U(N — 2) symmetry, while a normal solution 
lying on LC + and corresponding to u = preserves a £/(iV — 1) symmetry. For example, in 
the context of the two-Higgs-doublet model (N = 2) such a solution corresponds to a complete 
breaking of the electroweak symmetry group SU(2) x U(l). Such phase breaks the electric 
charge conservation, and makes the photon massive. 

Conditions when this phase appears were established in [28]. Since (r M ) that corresponds 
to such a non-symmetric phase is not restricted anymore to lie on the surface of LC + , the 
extremum condition of the potential takes a very simple form: 



(55) 



If B^y is non-singular, then solution of ( l55i) always exists and is unique. If the potential has 
any additional explicit symmetries, this symmetry is always conserved in this phase. If is 
singular, then depending on A^ Eq. (|55|) can have an empty set or a continuum of solutions. 

Whether the solution of (|55|) corresponds to a physically realizable extremum of the po- 
tential, depends on whether a M = {B~ l ) flv A v lies inside LC + . If it is so, then it can be a 
minimum or a saddle point. It is a minimum (and necessarily the global minimum) when B^ 
is positive-definite in the entire Minkowski space, i.e. when all Bi < 0. 

Search for the extrema on the forward lightcone LC + proceeds in the same way as before. 
One again introduces equipotential surfaces Ai c , but due to fixed eigenvalues Bi their geom- 
etry can be different. A typical MP can now be any 3-quadric: a 3-hyperboloid, a 3-ellipsoid, 
a 3-cone, or a 3-paraboloid. The geometric reformulation of the problem remains unchanged: 
the search for the global minimum corresponds to the search for the unique 3-quadric with 
the base point a M that touches but never intersects the forward lightcone. 

As a result, virtually all the statements about the number of extrema and minima, about 
the symmetries and their spontaneous violation remain the same. The only difference is that 
r M can shift from the surface of LC + inwards, and in order for an extremum on LC + to be 
a minimum, this shift must also increase the potential. It means that Lagrange multiplier A 
in ffl6|) must be positive. In fact, the eigenfrequencies uj 2 of oscillations that make iV-vectors 
ipi and ip2 non-parallel are proportional to A. In 2HDM, they correspond to the masses of 
charged Higgs bosons, [28 s [30] . 



10 Conclusions 

The aim of this paper was to provide an exhaustive description of the general two-order- 
parameter model with all possible ?7(l)-symmetric quadratic and quartic interaction terms in 
the mean-field approximation. The principal difficulty in the study of this model lies in the fact 
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that the Landau potential cannot be minimized with straightforward algebra. Here we showed 
that despite this difficulty, one can still learn a lot about the phase structure of this model. 
We developed the Minkowski-space formalism based on the reparametrization symmetry of 
the model and reformulated the minimization problem in simple geometric terms. We then 
proved several statements concerning the properties of the model (the number of extrema and 
minima, symmetries and their violation, the phase diagram). 

The most general 20P GL model can be viewed as a "template" for many particular re- 
alizations of the two-order-parameter model used in various condensed matter problems. We 
believe that by considering the most general case one can gain a more transparent understand- 
ing of phenomena taking place in particular situations, and one might even establish new links 
between seemingly unrelated models. 

We also note that the general method used in this paper (consider the model in the most 
general case, find the group of reparametrization symmetries, and using it find the structure 
behind the model) is very general and might turn out helpful in other circumstances. 

There remain several directions for future work. First, using dependence of the coefficients 
on temperature, pressure, etc., one can trace in detail the sequence of phase transitions as 
well as their critical properties in the mean field approximation. Second, one should study 
modifications caused by the presence of external fields (e.g. magnetic fields for two-gap super- 
conductors) and non-trivial boundary conditions. Third, one should analyze effects beyond 
the mean-field approximation, in particular, study how the symmetries of the model evolve 
under the renormalization group flow. Fourth, one should closely examine the existence, sta- 
bility and dynamics of the solitons. Finally, extension of the approach to models with several 
order parameters and/or with matric- valued OPs also appears to be feasible. 

I am thankful to Ilya Ginzburg, Otto Nachtmann and to the referees for discussions and 
useful comments. This work was supported by the Belgian Fund F.R.S.-FNRS via the contract 
of Charge de recherches and in part by grants RFBR 08-02-00334-a and NSh-1027.2008.2 



A Manipulation with 4-tensor B 



Here we collect some simple facts about the real symmetric 4-tensor B^ u . Let us first give 
explicit expressions for B^ u with raised or lowered indices: 
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(56) 



Here i,j = 1, 2, 3. Note that B^ v is not symmetric. 

Upon an 5*0(3) rotation, Boo remains invariant, while Boi and Bi d 
vector and symmetric 3-tensor, respectively. Upon a boost with "rapidity 
first axis, B^ transforms as: 
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where 



00 



«n = 



boo ~ hi . boo + hi , . , . , 

H cosh 2t] + boi sinn Zt] , 
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&00 + &11 
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cosh 2r/ + 6 i sinh 2r] , 



01 



sinh 2r] + b 01 cosh 2r] . 



0a 



b 0a cosh T) + b la sinh 77 , b' la = b la cosh 77 + b 0a sinh 77 , b' ah = b ab , a, 6 = 2, 3 , 



The eigenvalues -B, and eigenvectors of 5^ are defined according to 

B^e 1 ^ = B,i g^i) , or equivalently £/e w „ = B t e (i) „ . 



(58) 



The fact that .B^ is not symmetric implies that the some eigenvalues might be, in general, 
complex. However, as proved below, positive definiteness of B^ on the forward lightcone LC + 
ensures that they are real. 

In the diagonal basis one has: 
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If one considers a quadratic form in the space of 4- vectors p^ constructed on B^, then in the 
diagonal basis it looks as 

B^p» = B pl-J2B iP l 

i 

This quadratic form is positive definite in the entire space of non-zero vectors if and only 
if all Bi are negative. 



B Positive definiteness of V4 

A potential stable in a strong sense was defined as the one whose quartic part V4 is strictly 
positive definite in the entire space of the order parameters ipi except the origin. In the orbit 
space it corresponds to B^r^r" being positive definite on the entire forward lightcone LC + 
expect the apex. This criterion can be formulated in terms of the eigenvalues of B^ u : 

Proposition 4. Tensor B^ v is positive definite on the future lightcone expect the apex if and 
only if the following conditions are met: 

(1) B^ v is diagonalizable by an 5*0(1,3) transformation, 

(2) all spacelike eigenvalues Bi are smaller than the timelike eigenvalue B . 

Proof. Obviously, if B^ u satisfies conditions (1) and (2), then the positive definiteness follows 
immediately. So, one needs to prove that these conditions follow from the positive definiteness. 

The first step is to prove that the positive definiteness on LC + implies that all the eigen- 
values of B^ v are real. 
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Suppose, on the contrary, that there exists a pair of non-zero complex eigenvalues, b and 
b*, with respective (and necessarily complex) eigenvectors p M and g M : 

B» vV v = bp* 1 , B» v q v = &V . 

One can show that there can be only one pair of complex eigenvalues, thus, b is non-degenerate. 
Since B fl u is real, oc p^* (and can be taken equal to p^*). These eigenvectors are orthogonal, 
p^q^ = 0, which follows from the standard argument due to b ^ b*, and can be normalized so 
that p^Pn = q^q^ = 1. 

Consider now a non-zero real vector r^, 

= cp^ 1 + c*p* M , 

such that r'V^ = c 2 + c* 2 = 2|c| 2 cos(2</> c ) = 0. At fixed |c|, four such vector are possible. 
Take two of them: r^, 4 and r^ 4 , corresponding to C = 7r/4 and 37r/4. The quadratic form 
calculated on these vectors is 

B^r" = be 2 + b*c* 2 = 2\b\\c\ 2 cos (20 c + <f> b ) = T2\b\\c\ 2 sin(0 fe ) 

for ry A and r^ A , respectively. Since b is not purely real, sin(0 b ) ^ 0, in one of the two cases 
B^r^r" < 0, which contradicts the assumption. 

After all the eigenvalues of B^ v are proved to be real, the eigenvectors also can be chosen 
all real and orthonormal. These eigenvectors cannot lie on LC + (otherwise there would be a 
flat direction of V4), so there is one vector inside LC + with positive norm, norm p^po^ = 1, 
and three spacelike eigenvectors with negative norms p^Pi^ = —1 for each i = 1,2,3. Thus, 
the transformation matrix T that diagonalizes B^ v is real, and after diagonalization B^ u takes 
form diag(i?o, — B\, —B2, — -B3). Note that transformation T also conserves norm, so it can 
be realized as a transformation from the proper Lorentz group. 

Now, the requirement that B^ u is positive definite on LC + reads: 

Bq — (B\ sin 9 cos </> + B2 sin 9 sin <fr + B3 cos 9) > 
for all < 9 < n and <fi. This holds when Bq is larger than any Bi. □ 

Let us also see what changes for a potential stable in a weak sense. First, the statement 
that the eigenvalues are real and therefore eigenvectors can be also chosen real remains valid 
in this case. However, at least one eigenvector must now lie on the surface of LC + . This 
means that B^ v is in general not diagonalizable by the SO(l, 3) transformation group. More 
details are given in Section |U 



C Necessary condition for the spontaneous violation: 
explicit calculations 

Here, we show that the global minimum of the potential with all distinct Bi and = 
(A , Ax, A 2 , 0) can spontaneously break the Z 2 symmetry given by reflections of the third 
axis, only if B 3 is the largest spacelike eigenvalue: 

B 3 >B 1 ,B 2 . (59) 
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We assume, of course, that the vector lies inside the caustic cone: 



A 2 A 2 A 2 

+ 775 ^< ( p~ d~ \ 2 ( 6 °) 



(S 3 " ^O 2 (£? 3 " 5 2 ) 2 (So - #3 

This will be done by comparing the depth of the potential at the extrema that conserve and 
violate this symmetry. We will see that fl59l) is necessary for the pair of symmetry violating 
extrema to be the deepest ones. 

If (r M ) is an extremum point, then the potential at this point is 

V = -A,(r») + \B, u (r»)(r») = -\m^) = ~\b (r v ) ■ 



According to ffl8|) . the symmetry violating extrema take place at 
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The depth of the potential at this point is 
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Pick up another extremum (necessarily a symmetry-conserving one). It takes place at another 
r , which we rewrite as r = 7q • x. The depth of the potential at this point is 

\V\ = ^r (A - Ainx - A 2 n 2 ) 
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2 X 



A 2 A 2 A 2 

Ji a x /i 2 



B - B 3 {Bo - B 3 ) - (Bo - B x )x (B - B 3 ) - (B - B 2 )x 



Here, ni, n 2 are 



Til = ; r , Ho = ; : , + Ur, = 1 . 

A -{Bo-B x )ro' A - (B Q - B 2 )r ' 1 2 

Note that the last equation, in fact, is the fourth-order equation for tq. 

Difference between the two depths can be presented, after some algebra, in the following 
way: 

x)A ( A Aim A 2 n 2 



^ 3 ^ ^ 2 \ Bq — B 3 B\ — B 3 B 2 — B 3 ) ^ ^ 

The expression in brackets can be rewritten as a^n^, where n M = (1, n\, n 2 , 0) and 

a>=( A ° ^ A > ,0 

V B Q — B 3 Bx — B 3 B 2 — B 3 



From the caustic condition fl60l) one obtains a^a^ > 0, i.e. the four-vector a M lies strictly 
inside LC + . On the other hand, n M lies on the surface of LC + , and therefore, a^n^ > 0. 
Thus, the sign of the depth difference is given solely by the value of x. 
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If B 3 is the largest spacelike eigenvalue, then all symmetry-conserving extrema correspond 
to r"o < Tq , i.e. to x < 1. Therefore, all of them lie above the symmetry- violating points (and 
are saddle points, according to Proposition [2]). If B 3 is not the largest spacelike eigenvalue, 

(3) 

then there will necessarily be another extremum with ro > r$ , which corresponds to x > 1 
and, therefore, lies deeper that the symmetry-violating points (which are saddle points in this 
case) . 

An alternative, somewhat longer way to prove condition ( l59il using geometric properties 
of the potential was given in [29] . 



D Critical exponent for the distance squared function 
defined on a planar curve 

Let j(t) be a regular parametrization of a smooth plane curve, see e.g. [35]. Take a point on 
this curve, assuming that it corresponds to t = 0, and choose the coordinate frame at this 
point such that axis x is along the tangent and axis y is along the normal to the curve at this 
point, see Fig. [7J The curve then can be parametrized as j(t) = (X(t), Y(i)), with 

X(t)=t, Y(t) = 7 ^- + a n t n + o(t n ). (62) 

Here, Rq is the curvature radius of 7 at t — 0, while n > 2 describes the next higher order 
term. 




Figure 7: Regular curve j(t) and its evolute (dashed line); see text. 



Now, select a point r = (x, y) on this plane and calculate the distance squared from this 
point to the points of the curve, p 2 {t) = (X(t) — x) 2 + {Y{t) — y) 2 . This function has points 
of extrema at some values of t. For a generic point r, p 2 {t) will have a generic value at t — 0. 
However, if r*lies on the y axis, then p 2 (t) has a maximum or minimum at t = 0. At a special 
("critical") point along this axis, r — (0, Ro), 

p 2 (t) = X 2 (t) + (Y(t) - Ro) 2 « const - 2R a n t n + . 
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Therefore, the osculating circle has at 3-fold or 4-fold contact with the curve 7, depending on 
whether n = 3 or n > 4. As the point r moves along the y axis and passes through (0, Ro), a 
bifurcation takes place of the function p 2 (t). 

There are two possibilities to consider. If 7(0) is a generic point on a generic smooth curve, 
then expansion ( 1621) starts from n = 3, and at small t, p 2 {t) has a minimum and a maximum 
near at t = 0, both of which cannot be the global ones. Instead, if 7(0) is an apex of the 
curve (this is the situation shown in Fig. [7]), then expansion ( 1621) starts from n = 4, and at 
this critical point the minimum of p 2 {t) splits into minimum/maximum/minimum sequence. 
This is the only type of bifurcation the global minimum of the function p 2 (t) can participate, 
and we now focus on it. 

Let us shift f away from the critical point by a small amount, r = (e x , Rq + e y ), and 
recalculate p 2 {t): 

p 2 {t) « const + 1 4 ( 7^2 - 2a * R o J " - %- ■ (63) 

For a generic curve ^(t) (the second order curves included) the coefficient in front of t 4 is 
non-zero. Finding the minimum of (I63p . expanding p 2 near it, and extracting the coefficient 
in front of the quadratic term (t — t m i n ) 2 , which should behave as oc \e\ s , gives us the value of 
5. One can easily find that it depends on the direction of approach to the critical point: 

generic direction (e x 7^ 0) — > 5 = 2/3, 
symmetric approach {e x = 0) — > 8 = 1, (64) 

The latter case, by construction, effectively corresponds to the standard one-order-parameter 
Ginzburg-Landau model. 

This study can be generalized to the n + 1-dimensional case. Given a smooth n-manifold 
7(ti, . . . ,t n ), one can choose the coordinate frame a point of this manifold, that would corre- 
spond to the global minimum, align the coordinate axes with the eigenvectors of the quadratic 
term, and parametrize the manifold as 

X l = t l ,...,X n = t n , Y w + • • • + -|- + adi . (65) 

Here, we labeled the axes according to Ri < R 2 < ■ ■ ■ < R n , the bifurcation we study is at 
f = (0, . . . , 0, Ri). The calculations can be repeated giving the same result: if e x 7^ 0, then 
8 = 2/3, otherwise 8=1. Now e must be understood as the distance from the closest among 
the cusp points of the evolute. This result does not depend on the particular shape of the 
manifold, since it is essentially driven by the 4-th order nature of the bifurcation point. 
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